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ABSTRACT 

Space and astrophysical Alfven waves usually exist in a 
plasma containing thermal and nonthermal particles, and 
generally the thermal electrons and ions have unequal num- 
ber densities. Nevertheless, the linear properties of the 
waves are invariably assumed to be determined by a charge- 
neutral thermal plasma in the absence of the nonther- 
mal particles, while the nonthermal particles cause growth 
or additional damping superposed onto the background 
model. Here, both analytic plasma theory and numerical 
solutions of the dispersion equation are used to show that 
this prescription is seriously flawed even for stable plas- 
mas; even when the nonthermal particles do not contribute 
significantly to the dispersion equation, unequal thermal 
electron and ion number densities (due to the presence 
of the nonthermal particles) may cause fundamental low 
wavenumber modifications to the Alfven modes, including 
the creation of a new resonance and severely modified dis- 
persion. These results are found for both cold and warm 
plasmas. Previous work on Alfven waves should be re- 
evaluated in view of these results. 

Keywords: Alfven Waves, Dispersion Relations, Kinetic 
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1. INTRODUCTION 

Alfven and fast magnetoacoutic waves are the basic low 
frequency (u; < 0*) electromagnetic waves in a magnetized 
plasma. These waves are of great importance in many fields 
of space physics and plasma astrophysics, including parti- 
cle acceleration at solar system and astrophysical shocks 
(Ref. 1), the production and maintenance of cosmic rays 
(Ref. 2), turbulence theory and nonlinear waves such as 
solitons. Recent research has generally been concerned with 
the nonlinear properties of, and detailed applications in- 
volving, Alfven waves rather than on the linear properties 
of the waves such as mode structure and dispersion relations 
which have been assumed well established. In contrast the 
linear properties of Alfven waves are explicitly considered 
in this paper, and important new results are found which 
imply a need for re-assessment of much work on detailed 
applications and nonlinear physics involving Alfven waves. 

Alfven waves in space and astrophysical applications typi- 


cally exist in plasmas containing both thermal and nonther- 
mal particles, and generally the thermal electrons and ions 
will have unequal number densities due to a preponderance 
of one energetic species. The linear properties of the waves 
Eire, however, generally described using the following “con- 
ventional” prescription: dispersion properties (such as the 
mode structure and dispersion relations) are determined in 
a charge-neutral thermal plasma in the absence of the non- 
thermal particles, while the nonthermal particles merely 
cause growth or additional damping superposed onto the 
background plasma modes. This prescription is shown to 
be seriously flawed at low wavenumbers in this paper. Mod- 
ifications to the linear properties of the waves, due to the 
presence of nonthermal particles are found to arise in two 
ways: (1) Unequal thermal electron and ion number den- 
sities cause fundamental low wavenumber modifications in- 
cluding the creation of a new resonance and severely modi- 
fied dispersion (Ref. 3). (2) Direct contributions to the wave 
dispersion equation by the energetic particles may modify 
the Alfven waves at low wavenumbers (as implied by the 
resonance condition u zb £l a — k\\V\\ where a denotes the 
particle species). 

In this paper I consider the effects of unequal thermal num- 
ber densities in some detail, using firstly analytic plasma 
theory and then numerical solutions of the full wave dis- 
persion equation. The direct effects of the energetic par- 
ticles on the wave dispersion will be considered elsewhere 
(e.g., Ref, 3); note that others (Refs. 4, 5) have previously 
considered such effects. 

2. MODIFICATIONS DUE TO UNEQUAL THERMAL 
NUMBER DENSITIES 

The dispersion equation for parallel propagating electro- 
magnetic waves in a magnetized plasma with both nonrel- 
ativistic thermal (Maxwellian, T) and 

k 2 ~ w 2 ~Y,%- = 

T 1 NT * 

with £ a = W/y/2KV a , (Z = (W + <7 e a n a /n i )/V2KV a , 
V a is the thermal speed divided by c, and Z((%) is the 
Friede-Conte function. Here W = K = kc/Qi , a 

is the handedness (±1 or RH/LH) of the mode, and e ay 
il ay u> pay and F a (W y K) are the charge (unsigned) gyrofre- 
quency, plasma frequency and plasma dispersion function, 
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respectively of species a. Wave frequencies well above rele- 
vant collision frequencies have been assumed. For nonrela- 
tivistic particles with isotropic distribution function / a (p), 
F q (W,K) is proportional to 



p\ l-t* 2 ) 

(W + ae a Sl a /n t -Kpp) 


dfa 

dp' 


( 2 ) 


where p is the particle momentum normalized by c and p 
is the cosine of the pitch angle. Accordingly, each F a de- 
pends on the distribution f a of nonthermal particles and 
may differ considerably from the thermal plasma contribu- 
tion to the first sum in (1) for the same number density of 
particles. 

In cold or warm thermal plasmas where V a < 10 -2 , the 
large argument expansion of the Friede- Conte function is 
rigorously justified for Alfven waves with If < 1 when the 
Alfven speed I £ is much less than c, leading to the so-called 
“cold plasma” description of the terms in the thermal (T) 
sum in (1): 




w 


ft? W + <re a n a /Qi 




NT 


(3) 


An analytic description of the low wavenumber modifica- 
tions due to unequal thermal electron and ion number den- 
sities, and the circumstances under which these effects may 
be relevant, proceeds as follows. Consider a charge- neutral 
electron-proton plasma with cold or warm (V Q < 10" 2 ) 
thermal components and a single (for simplicity) nonther- 
mal component, with total ion and electron number densi- 
ties n 0 - After minor algebra and the subsequent assumption 
W <C 1 (w in the denominator of the thermal sum 

(T), (3) becomes 


K 2 -W 2 - 


m e 


m;ft 2 
- Fnt — 0. 


W 2 (u 2 pe +u> 2 pi ) + aW(« 



(4) 


Conventional derivations assume that the thermal elec- 
tron and ion number densities are equal, so that the W 
term in the square brackets vanishes, and that the non- 
thermal term may be ignored. Using the usual definition 
of the Alfven speed V A {c 2 /Vf = 311(1 writ ' 

ing «£.,,• = n 0 e 2 /m,< 0 , one finds that c 2 /\tf = 1 /U? = 
u po ,-(l + n»«/mj)/ft? and the conventional derivation yields 

K 2 - W 2 - W 2 /U a = 0, (5) 


which for U A < 1 (i.e., V A < c) gives W = KU A , or w = kV A 
as required. 

Here, the electron and proton thermal number densities are 
permitted to differ from each other, as expected in a typical 
space or astrophysical plasma, so that the W term in the 
square brackets in (4) is retained: this term is responsible 
for the low wavenumber modifications to the Alfven modes 
now described. 

The effects of unequal thermal particle number densities on 
the wave dispersion are clearly separated from the effects of 
direct contributions (through F NT ) to the dispersion equa- 
tion by the nonthermal particles when the direct contribu- 
tion Fnt of the nonthermal particles is negligible, i.e., the 


form of the distribution /att is such that F is much less than 
the “cold plasma” contribution of the corresponding num- 
ber density of thermal particles removed from the thermal 
sum in (4). Then, assuming that 6n = n e — n; (the dif- 
ference in thermal number density) is small compared with 
no and that ft < 1, the dispersion equation (4) becomes 

W 2 + a — W - K 2 ltf = 0. (6) 

no 


Below a distribution function f Q suitable for describing non- 
thermal particles accelerated at shocks is shown to satisfy 
this condition that Fnt be ignorable. 

For real W (i.e., neglecting damping) the solutions of (6) 
are 



The analytic predictions therefore are: (1) At high wave- 
numbers K such that KUa > |£n|/no, i.e., the third term 
in (6) dominates the second term, the ordinary Alfven 
dispersion W ~ KUa is recovered. (2) Fundamental low 
wavenumber modifications occur when KUa < \Sn\/riQ. In 
particular, (3) the mode with a6n < 0 has a new resonance 
at W ~ \6n\/n 0 . (4) The other mode approaches a disper- 
sion relation W ex K 2 . In summary, the effects of unequal 
thermal particle number densities may cause fundamental 
low wavenumber modifications to the Alfven modes; the di- 
rect contributions Fnt of nonthermal particles may either 
aid or compete with these modifications due to unequal 
thermal particle number densities. 

In space and astrophysical plasmas shocks axe observed 
(Ref. 1) or inferred (Ref. 2) to accelerate nonthermal parti- 
cles with distribution functions f(p) oc p~ a with a ~ 4. It 
is therefore natural to consider nonthermal particles with 
distribution function 

/(p) “ (p 2 + q 2 ) 2 (8) 

for p < p m ax and zero elsewhere, in this connection, 

since this distribution is analytic with /(p) oc p~ 4 at high 
momenta p > q. The plasma dispersion function Fnt for 
this distribution has been calculated analytically using the 
Plemelj formula in (2) [justified when the imaginary part of 
W is much less than the real part, as found in the following 
examples]. Numerical solutions to the dispersion equation 
(3) for a three- component electron- proton plasma compris- 
ing cold electrons and ions and nonthermal particles (with 
distribution function (8) for p max =0.1 and q = 10“ 2 ) are 
now presented. 

Figure 1 shows the W-K dispersion curves for the above 
charge- neutral plasma with = 300 ~ c/\J contain- 

ing various relative number densities of nonthermal elec- 
trons, i.e., deficits of thermal electrons 6n < 0. The pure 
cold plasma results are also shown for comparison. The di- 
rect contribution Fnt(W, K) of the nonthermal particles to 
(3) is found to be numerically ignorable, thereby fulfilling 
the requirements of the analytic theory above. Both the low 
and high wavenumber regimes are shown for the <r = -fl 
(full lines) and a = - 1 (dashed) modes. At high wavenum- 
bers the curves closely follow the pure cold plasma solutions 
as the RH fast magnetoacoustic (<7 = +1) and LH Alfven 
(cr = —1) modes become the whistler and ion cyclotron 
modes, respectively. Deviations from the pure cold plasma 
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curves become apparent at wavenumbers KVa ~ |6n|/n 0 . 
Once K < \Sn\/noUA , the a = +1 mode approaches a reso- 
nance at W ~ |6n|/n 0 , while the a = — 1 mode approaches 
a dispersion relation W oc K 2 . The situation for nonther- 
mal ions is shown in Figure 2; it is identical to Figure 1 
except "that the handedness of the modes approaching the 
new resonance and dispersion W oc K 2 are reversed, con- 


A-G8B-SI3 


NON -THERMAL ELECTRONS 



Figure 1. Dispersion curves for a cold plasma with various 
nonthermal electron densities for RH (full lines) 
and LH (dashed lines) waves verifying the ana- 
lytic predictions. At low wavenumbers the pure 
cold plasma solutions are denoted by the lowest 
full line. 


A - 608 - 514 



Figure 2. As for Figure 1, except for nonthermal ions. 


sistent with (7). Thus, Figures 1 and 2 verify all the predic- 
tions of the analytic theory derived above, in a physically 
relevant setting. In addition, the full dispersion equation 
(1) has been solved for the nonthermal distribution (8) for 
a wide range of plasma temperatures: for V a < 10 -2 the 
cold plasma description is found to be accurate to within 
1 part in 10 8 for frequencies W < 0.1. Thus, the effects of 
unequal thermal electron and ion temperatures described 
above should exist and be potentially important in plasmas 
with either cold or warm thermal components. 

3. SUMMARY AND FURTHER WORK 

In summary, using analytic kinetic theory and numerical 
solutions of the wave dispersion equation, I have shown 
that: (1) Unequal thermal electron and ion number densities 
(due to the presence of nonthermal particles) in a charge- 
neutral plasma may lead to fundamental low wavenumber 
( K < \Sn\/n 0 UA) modifications of the Alfven modes, in- 
cluding (2) the creation of a new resonance at W ~ |£n|/n 0 
for the mode with aSn < 0 and dispersion W oc K 2 for the 
other mode. The nonthermal particles may also directly 
modify the Alfven modes (e.g., Ref. 3) through the terms 
F a in (1); these direct modifications aid or compete with 
the effects of unequal thermal number densities. The ba- 
sic requirement for the effects of unequal thermal number 
densities to be significant in determining the wave disper- 
sion is that the direct contributions F a of the nonthermal 
particles be small (due to the form of the particle distri- 
butions f Q ) compared to the “cold plasma” contribution of 
an equal number density of thermal particles. This require- 
ment should not be considered too stringent: firstly, it is the 
conventional approximation made in basically all previous 
discussions of space and as trophy si cal Alfven waves, and 
secondly, a natural, analytic distribution mimicing the 
(nonrelativistic) power-law distributions expected from 
shock acceleration in space and astrophysical plasmas obeys 
this requirement. Thus, in view of the expected ubiq- 
uity of physically interesting plasmas containing unequal 
thermal electron and ion densities, much previous work 
on detailed applications and nonlinearities involving Alfven 
waves should be re-evaluated in view of these results. 

One immediate implication for studies of Alfven wave non- 
linearities may be seen in Figures 1 and 2: conventionally 
Alfven waves are considered modulationally unstable where 
d 2 u>/dk 2 < 0, as for the LH (<j — — 1) mode but not the RH 
(<r = -hi) mode in a pure cold plasma at high wavenum- 
bers. Figures 1 and 2 indicate that the new low wavenum- 
ber resonance, due to the effects of unequal thermal particle 
number densities, is modulationally unstable according to 
this criterion. Moreover, in contrast to the conventional 
picture, these low wavenumber resonances are undamped 
even in warm thermal plasmas, and either the LH or RH 
mode may be modulationally unstable depending on the 
species with a deficit in thermal number density. Further 
research on this matter is needed, a s well as on the detailed 
connection between this kinetic theory and MHD theory, 
the effects of additional ion species, and other distributions 
of nonthermal particles. 
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